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In this work we extend our previously developed formalism of Newtonian multi-fluid hydrody- 
namics to allow for coupling between the fluids and the electromagnetic and gravitational field. 
This is achieved within the convective variational principle by using a standard minimal coupling 
■ prescription. In addition to the conservation of total energy and momentum, we derive the conser- 

vation of canonical vorticity and helicity, which generalize the corresponding conserved quantities 
^vq , of uncharged fluids. We discuss the application of this formalism to electrically conducting systems, 

magnetohydrodynamics and superconductivity. The equations of electric conductors derived here 
are more general than those found in the standard description of such systems, in which the effect 
of entrainment is overlooked, despite the fact that it will generally be present in any conducting 
multi-constituent system. 
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Ch ' I. INTRODUCTION 

>>: 

In the previous paper pj, henceforth referred to as Paper I, we have shown how the convective variational principle 
can be used to derive the general equations of motion for a system of interacting uncharged fluids. This variational 
principle is also at the heart of a series of papers by Carter and Chamel based on a fully covariant spacetime 

! formulation of Newtonian hydrodynamics, which is formally closer to the relativistic formulation. 

Here we extend our "3+1" framework developed in Paper I to allow for charged fluids and their coupling to the elec- 
tromagnetic and gravitational field. In Paper I we already included the gravitational field via an explicit prescription 
£>v for the external force, but now this coupling is derived in a more natural way from the variational principle using the 
' same minimal coupling procedure as for the electromagnetic field. 

Due to the coupling to the electromagnetic field, the resulting theory can only be considered an "approximate" 
Newtonian framework, as strict Galilean invariance will be violated. While the Newtonian hydrodynamic equations 
t-H , are strictly invariant under Galilean transformations, the equations of electrodynamics are invariant only under 
Lorentz-transformations. This well known discrepancy lead to the development of relativity, and a coupling between 
Newtonian physics and electromagnetism is strictly speaking either inconsistent or selects a preferred frame ("ether") . 
As shown in , demanding strict Galilean invariance of electrodynamics without an ether forces one to adopt either 
an "electric" or "magnetic" limit of the theory, in which certain essential effects would be absent (e.g. the magnetic 
force between electric currents or local charge conservation). The only fully consistent approach is to work within a 
I/"} (locally) Lorentz-invariant relativistic framework, as used for example in [if 0- 

From a practical point of view, however, a "3+1" non-relativistic formalism is often fully sufficient in terms of 
precision and usually more easily applicable to many problems. The aim of this paper is therefore to provide a 
flexible and general framework for the description of a wide variety of charged multi-constituent systems in the "non- 
relativistic" regime of small velocities and low frequencies (c -2 ) 1 . We emphasize that we do not attempt to construct 
a strictly Galilean- invariant theory, rather this should be regarded as a truncated theory at order (c~ 2 ) of a Taylor 
*-£h . expansion of the underlying fully covariant theory. 

For simplicity we restrict our analysis to non-magnetic and non-polarizable fluids, such that the vacuum Maxwell- 
equations retain their form on the macroscopic level of hydrodynamics, and the interaction between matter and fields 
is restricted to the minimal coupling type. The inclusion of electric and magnetic polarization is postponed to future 
work. 

Using this variational framework, the fundamental effects of "entrainment" and "chemical coupling", inherent to 
multi-constituent systems (as discussed in Paper I), are automatically included in the formalism. This is a substantial 
improvement over the standard "orthodox" description of electrical conductivity and charged fluids in general, in 
which these effects are usually completely overlooked. The present framework is also more general in the sense of 
being readily applicable to an arbitrary number of interacting constituents and fluids. 
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The plan of this paper is as follows: in Sec. II we extend the variational principle introduced in Paper I to the 
case of ("separable") coupling to the electromagnetic and gravitational fields, and derive the equations of motion for 
this coupled system. In Sec. Ill we derive the conservation of charge and mass as well as of energy and momentum 
from the variational principle. Sec. IV is devoted to conserved quantities under transport of the fluid flow, namely 
the canonical vorticity and helicity, and their special relation to superconductors. In Sec. V we discuss applications of 
this formalism to particular physical systems. In the appendix we derive the non-relativistic form of electrodynamic 
equations and show their "approximate" Galilean invariance. 

II. VARIATIONAL DESCRIPTION OF MULTI-CONSTITUENT SYSTEMS 

A. Kinematics 

We follow the notation introduced in Paper I, the various constituents are indexed by capital letters X, Y, which 
run over the constituents labels. The fundamental kinematic variables are the constituent number-densities nx and 
the associated transport currents nx, which are related to the respective flow- velocities vx as nx = nxVx- As in 
Paper I, we define the relative velocities Ajy between fluids X and Y as 

Axy = v x - v Y ■ (1) 

In addition to the masses per particle m x , we now admit another transported quantity, namely the charge per particle 
q x . The corresponding total densities and currents are therefore the mass density p and mass current p, given by 

p = ~^^nxm x , and p = ^m x n^, (2) 

and also the charge density a and charge current j, defined as 

g = ^ n x q X , and j = ^ n x q X , (3) 

where here and in the following ls used to indicate the sum over repeated constituent indices, no automatic 
summation over constituent indices is assumed. 

In addition to the dynamics of the matter-system we also want to include the coupling to the gravitational gauge 
field $ and the electromagnetic gauge fields Ao and A. The corresponding gauge-invariant field quantities are the 
gravitational acceleration g, which is 

fl = -V$, (4) 
and the electric and magnetic fields E and B, which are 

E = VA Q - -d t A , and B = V x A . (5) 

c 

These fields are invariant under the gauge transformations 

*' = * + C(t), A ' = A + -d t i> , and A' = A + Vip. (6) 

B. Dynamics 

The dynamics of the system is described by an action X of the form 

1 = J KdVdt, (7) 

in terms of the (generally gauge-dependent) total Lagrangian A. Using the standard minimal coupling prescription, 
and assuming the electrodynamic field-contribution to be fully separable from the hydrodynamic Lagrangian, we can 
postulate the form of A to be 

A = A H {n x , n x ) + A EM (E, B) + A grav (g) + (aA + -J ■ A) - p<$> . (8) 
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The gravitational field Lagrangian A grav has the well-known expression 

W = -g^(V*) a , (9) 

where G is Newton's gravitational constant. The electromagnetic field contribution (cf. appendix) in the non- 
polarizable case is given by 

A EM = i- (E 2 - B 2 ) . (10) 

As we have seen in the case of uncharged fluids (cf. Paper I) , the hydrodynamic Lagrangian A H defines the dynamical 
momenta p x and p x per fluid particle as 

dA H = J2(Podn x +p X >dn x ) , so p x = |^ , p x = J^2- . (11) 

In a similar manner the variation of the electrodynamic field Lagrangian A EM defines the so-called "electric displacement 
field" D and the "magnetic field strength" H as conjugate variables to the electromagnetic fields, namely 

dA EM = ^-D ■ dE - ^-H ■ dB , (12) 
47r Ait 

and using the explicit form iflfljl we find 

,-, . <9A EM dA EM 

D = 4n— = E, H= 4n— = B. (13) 

Although we can trivially identify D = E and H = B in the present non-polarizable case, in the following we 
nevertheless keep the formal and conceptual separation between the "kinematical" E and B and the "dynamical" (i.e. 
derived from A) quantities D and H. 

Due to the presence of the minimal coupling terms in © , the variation of the total Lagrangian A generalizes the 
dynamical momenta p x and p x to the (gauge-dependent) canonical momenta tt x and 7r x , namely 

dA = ^2 ( w o dn x + ^ X ■ dn x ) _ P d ® + G dA o + ~3 ■ dA + rfA EM + dA grav , (14) 
from which we obtain the relations 

*o = Po +q X A ~m x <i>, (15) 

n x = p x + l q x Aj n 6 ) 

c 

expressing the canonical momenta in terms of the dynamical momenta and the gauge fields. 

C. The equations of motion 

We see from definition (0 that the first two Maxwell equations are satisfied identically, i.e. 

V • B = , and V x E + -d t B = . (17) 

c 

Using 11211 and H14I1 it is not difficult to show that requiring invariance of the action under free variations 5Aq and 
5 A of the gauge fields results in the two remaining Maxwell equations 

1 47T 

V • D = 47TCT , and V x H - -d t D = —j . (18) 

c c 

Similarly, requiring invariance of the action with respect to free infinitesimal variations <5$ straightforwardly results 
in the Poisson equation for the gravitational field, 



V 2 $=47tG/9. (19) 
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The derivation of the equations of motion for the fluids is formally completely equivalent to the uncoupled case 
discussed in Paper I. The only difference is that the dynamical momenta p x and p x are now replaced by the 
corresponding canonical momenta w x and tv x . 

One has to consider infinitesimal spatial displacements £ x and time-shifts tx of the flowlines, which should be 
considered the "true" fundamental quantities describing hydrodynamics (corresponding to the Lagrangian framework) , 
while free variations for densities and velocities (characteristic for the Eulerian framework) fail to produce the correct 
equations, except when adding "ad-hoc" constraints to the Lagrangian. The Eulerian variables, density nx and velocity 
vx can be expressed in terms of the underlying flowlines (and initial conditions) alone, and one should therefore also 
consider the hydrodynamic Lagrangian A H in © as a function of these variables. These subtleties relating to the 
hydrodynamic variational principle are discussed in greater detail in Paper I. 

The resulting induced variations of the densities, Snx, and currents, Snx have been derived in Paper I as 

5n x = -V • [n x £ x ] + \ n x • Vt x - T X d t n x ] , (20) 
Sn x = nxdtZx + (nx-V)Zx-(Zx-V)nx-nx{V-Zx)-dt(n X Tx) ■ (21) 

Substituting these expressions into the variation of the total Lagrangian 11411 and integrating by parts, we arrive at 
the following form of the induced action variation 



SI 



jY.^ X rx-f X -ix)dVdt, (22) 



where the canonical force densities f x (acting on the constituent) and the energy transfer rates g x (into the con- 
stituent) are found explicitly as 

f X = n x (d t iv x - Vtt x ) - n x x (V x tv x ) + n x r x , (23) 

g x = v x ■ (f x - n x r x ) - 7T x r x , (24) 

where rx is the particle creation rate for the constituent X, i.e. 

r x = d t n x +V-n x . (25) 

The canonical force density f x is the total (canonical) momentum change rate of the constituent X , and the last term 
in l(23jl represents a contribution that is purely due to the change of the particle number. By inserting the explicit 
expressions Ill5llfijl of the canonical momenta into l(23"jl and il24"jl . we can separate out the purely hydrodynamic 
contribution f x , which expresses fluid inertia and pressure, and which has the same form as in the uncharged case 
(Paper I), namely 

f x =n x (d tP x -Vp x ) -n x x {V xp x ). (26) 
Introducing the gravitational and electromagnetic forces f x leiv and f x M , defined by their usual expressions 

ifrav = ~n x m x V^ 1 (27) 
/™ = n x q X (e + -vx x b] , (28) 



we can now rewrite the canonical force densities l|23j) in the form 

f X = fn - Z^-av - /em + r X !T X . (29) 

The actual physical equations of motion are obtained by prescribing the canonical force densities f x acting on the 
fluids. If we only require invariance of the action for a common displacement and time shift £ x = € an d tx — t, we 
obtain the minimal equations of motion for the total system, namely 

^/ X = / cx t, and Y,9 X =9c*t, (30) 

where / ex t and <7 ex t are interpretable as the "external" force density and energy transfer rate acting on the system. 
This generalizes the more common action principle of isolated systems, in which the external influences / ex t and g ox t 
vanish and therefore the equations of motion would be obtained by requiring the action to be invariant under small 
variations. "External" here is meant in the sense of not being included in the total Lagrangian, which could also 
include, for example, viscous forces. 
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D. The hydrodynamic Lagrangian A H 

As shown in Paper I, the hydrodynamic Lagrangian density A H for the class of "perfect" (i.e. with an isotropic 
energy function £) multi-fluid systems is given by 

A H ( nx ,%) = Vm^-^x,A^), (31) 

where A^y denotes the relative velocity between fluids X and Y as defined in QJ. The total differential of the 
thermodynamic potential £(nx , A^ y ) determines the first law of thermodynamics for the matter-system, namely 

d£ = ^2fi x dn x + - aXY dA XY > (32) 

X.Y 

which defines the chemical potentials /j, x and the symmetric entrainment matrix a XY . The dynamical momenta p x 
and p x defined in l|TT)) are therefore found as 

x x 2a XY , . 

p x = m x v x - > Axy , 33 

-Po = fi x -m x ^-+v x -p x . (34) 

III. GAUSS-TYPE CONSERVATION LAWS 

A. Conservation of charge and mass 

We see from 1)26)1 - 1128)1 that the force contributions f x , f x rav and f x u are invariant with respect to gauge trans- 
formations of the gravitational and electromagnetic fields. However, the last term in the expressions 1(23)1 and l(24*)l 
(accounting for momentum and energy change due to particle number changes) is generally gauge dependent. While 
this is not prohibited for individual constituent forces, the total equations of motion 1)3011 have to be gauge invariant, 
and by using II15I16I) we can therefore deduce the constraints 

Y J 1 Xr x = d t <j + V-j = 0, (35) 
Ym x r x = d t p + V-p = 0, (36) 

where the total densities and currents have been defined in <|2j and . Gauge invariance therefore implies conservation 
of the "charge" associated with the gauge field. As a further consequence we find the following useful relations, 

J2r X 7T X =Y J TxP X , and J2r x n x =J2 r xPo ■ (37) 

We note that the conservation of mass was already derived in Paper I as a consequence of Galilean invariance, which 
can also be considered a gauge freedom. 

B. Momentum conservation 

As shown in Paper I, the purely hydrodynamic force densities f x satisfy the relation 

1 + r x P x l ) = d t x n + VjTi? , (38) 

where the hydrodynamic momentum density J H and stress tensor are given by 

J H = 5> x p X , and IV »'y/' V ' ' *.'/"• (39) 
The "generalized pressure" is defined via the Legendre transformation of A H , namely 

<J = A H - ^ ( n xPo + nx ■ P X ) , (40) 
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and gij are the components of the metric tensor, which in Cartesian coordinates is simply <7y = 5ij . Using these 
relations together with ll29)) and QJJ, we can write the total force balance equation ffifijl in the form 

i^xt-E^^^+W -^ra V -/^M, (41) 

where / gra v = fa&v 

tillQ J" EM = J2 /em are the total gravitational and electromagnetic force densities. Using Q3J 
and 1(281) we find explicitly 

/grav - ~pV$, (42) 

/em = crB + x B . (43) 
One can easily verify that the gravitational force term can be written as the divergence of a tensor, namely 

- /* rav = y 3 T^ , with T« av - ^-L f V W<I> - i (V$) 2 . (44) 

Using the Maxwell equations II17J) , II18II and the total differential 11211 of A EM , we can show that the total electromag- 
netic force density lH3)l can be similarly rewritten as 

-fL=dtJ*J+V j TH i , (45) 

in terms of the momentum density J EM of the electromagnetic field, 

Jem = 7— -D X B , (46) 

and the Maxwell stress tensor T^m, which is found as 

= -i- (£W + B^) + (a em + ^B • b) fl « . (47) 

Putting all the pieces together, we obtain the following form for the total momentum conservation 14111 : 

d t (J* + J EM 4 ) + V,r^ = /* xt , (48) 

where the total stress tensor is given by 

rpij _ ™ , /•,./ , ™ 1 !(); 
L H ' ^grav ^ 2 em ■ \^^J 

An important property of the total stress tensor T 13 is that it is symmetric. The symmetry of the gravitational part 
(j4*4l is obvious, while the symmetry of the hydrodynamic stress tensor ijSfljl has been shown in Paper I. It remains to 
prove the symmetry of the electromagnetic stress tensor 114 7J1 . In the present non-polarizable case, this follows trivially 
from the identities Ill3|) , It is interesting to note, however, that this can also be derived more generally as a Noether 
identity of the variational principle, assuming only the separable form (JHJ of the Lagrangian. 

In order to show this, we extend the variational principle slightly by admitting also metric variations Sgij, so that 
Ill2t now reads as 



<5A EM = 7— B • SE - y-B ■ SB + 6 gij . (50) 

4-7T 47r ogi 



Consider an active time-independent infinitesimal displacement £ of the whole system including the metric, which 
induces the following Lagrangian changes: 

AEi=E l V&, AB t = B l V&, Afly = -2V (i £j) . (51) 

Using these transformations together with J 5 Oil , we obtain the induced Lagrangian change of A EM as 

9A EM ~ 



1 

AA EM = — 

47T 



D l E 3 - B l B 3 - 8tt- 



. (52) 
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This active transformation is equivalent to a coordinate transformation — £, and therefore the requirement of A E 
being a scalar is AA EM = 0, which leads to the associated Noether identity 



2 



<9A EM 



-3- [D* E 1 - H l B 1 ] . (53) 



dg lj An 

Using the manifest symmetry of the left-hand side, we obtain 

E l D j + //' B j = E ] D i + II ! IV , (54) 

which concludes the proof of the symmetry of T EM . Note that in general Tem need not be symmetric, only the sum of 
all contributions to T %3 is subject to this constraint. The symmetry of T E „ is a special consequence of the assumption 
of a "separable" interaction of the form © . 

C. Energy conservation 

We have seen in Paper I that we can write 

{vx ■ Si - r x p*) = d t E H + V • Q H , (55) 

in terms of the hydrodynamic energy density E H and energy flux Q H , which are given by 

E H = Y j nx -P X - Ah, and Q H = ^(-p$)n x , ( 56 ) 

while for the gravitational and electromagnetic work contributions <? gr av = v x • /^ a v an d 9em = Tl v x ■ /f ul we 
find using l(57jl and l(28 ]l : 

ffgrav = -p-V$, (57) 

5em = j ■ E . (58) 

Using these expressions, the energy equation l(3"0j) can be written as 

5cxt = d t E n + V • Q H - ffgrav - ^em ■ (59) 

Using Maxwell's equations l|17ll and l|18|l . one can write the electric work g EM in the form of a conservation law, namely 

- 5em = d t E EM + V • S , (60) 
where the electromagnetic field energy density E EM is given by 

E EM = ^—E ■ D A EM = — -— (E 2 + B 2 ) , (61) 

47T 07T 

and the second equality was obtained using the explicit Lagrangian 11011 . The energy flux S is given by the Poynting 
vector 

S=—ExH. (62) 

In the present non-polarizable case, i.e. D = E and B = H, we recover the well-known relation between the energy 
flux and momentum of the electromagnetic field, namely J EM = S/c 2 . Summarizing, we can cast H59H in the form of 
a conservation of total energy, namely 

d t (E n + E EM ) + V • (Q H + S) = 5cx t - p ■ V$ . (63) 

We note that formally one can also write the gravitational work in the form of a conservation law, but the expression 
for energy density and flux are neither unique nor gauge invariant, and one can also not eliminate the mass current 
p from these expressions due to the lack of a dynamic law for the gravitational field in the Newtonian framework. 
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IV. CONSERVATION ALONG FLOWLINES 



In this section we show how the conservation of vorticity and helicity, derived for uncharged fluids in Paper I, can 
be generalized quite naturally to the case of fluids coupled to the electromagnetic and gravitational field. We note 
that the technical steps involved in this discussion are largely analogous to the treatment in Paper I, and we therefore 
skip most intermediate steps. 



A. Generalized Kelvin-Helmholtz vorticity conservation 

We define the hydrodynamic vorticity 2-form w as the exterior derivative of the dynamical momentum 1-form p, 
namely w = dp, and the more common dual vorticity vector W, which is W = V x p. 

In the presence of electromagnetic fields, the more fundamental quantity is the canonical vorticity 2-form w_, which 
is defined in the same way but with respect to the canonical momentum n, namely 

w = dn_ , (64) 

and the dual canonical vorticity vector W is therefore given by 

W = V x 7T . (65) 

With we see that the relation between canonical and hydrodynamic vorticity is simply 

w = w+-dA, and W = W + - B . (66) 
c c 

We note that by the Poincare property (namely dd — 0) , the exterior derivatives of the vorticity 2-forms vanish 
identically, i.e. dm = 0, which equivalently expresses the fact that the vorticity vectors are divergence-free, i.e. 

vw = o. 

We can write the expression i'2'Mi for the canonical force / acting on one constituent in the language of forms as 

d t TL + v\ dw - dn = - (/ - Ttt) , (67) 
n — 

where J indicates summation over adjacent vector- and form-indices. 

In the following it will be convenient to separate the "proper force" per particle acting on the right-hand side of 
H67H into its non-conservative part 3; and a conservative contribution d<j), namely 

-(/-Ar) = + (68) 
n — ~ 

Applying the Cartan formula for the Lie derivative of a p-form to the 1-form tv, namely £ v n = v\ dn_ + d(v\ n) , allows 
us now to rewrite the force equation l|67H more conveniently as 

(d t + £ v )Tr = dQ + $ 1 (69) 

where the scalar Q is given by Q = ttq + v\n + <f>. Lie derivatives and partial time derivatives commute with exterior 
derivatives, so we can apply an exterior derivative to l(69"jl and with l(64Tl we obtain the Helmholtz equation of transport 
of canonical vorticity, namely 

{d t + £ v )m = d$, (70) 

which shows that the canonical vorticity is conserved under transport by the fluid, if the proper force per particle 
acting on the fluid is purely conservative, i.e. if # = 0. In its more common dual form, this equation can be written 

as 

9 f W-Vx(«xW) = VxJ, (71) 

Substituting the explicit relation l(66)l between the dynamical and the canonical vorticity, and using Maxwell's equation, 
this can be re-expressed as 



d t W -Vx{vxW) = Vx 



J? + q [E+-VXB 



(72) 
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which shows that in the case of a charged constituent, the dynamical vorticity W is generally not conserved even in the 
absence of a non-conservative external force due to the presence of the Lorentz force acting on the fluid. However, 
the canonical vorticity is conserved in this case and therefore generalizes the vorticity conservation of uncharged fluids. 
For the canonical circulation C of a closed circuit <9E, which is the boundary of a 2-surface E, we have 



C= f tt = / w = / W ■ dS , (73) 

J <9£ JT. JT, 

where dS is the surface normal element. We see from ijfifif that the canonical circulation C can also be expressed as 
the sum of the hydrodynamic vorticity flux (i.e. dynamical circulation) and the magnetic flux through the surface E, 
namely 

C= t W-dS+- [ B-dS. (74) 
For the comoving time derivative of the circulation C we find using ij69l) 

Tt = f ( " dt+£v) - = i^- (75) 

which is Kelvin's theorem for the conservation of canonical circulation. We note that strict conservation only applies 
if the non-conservative force per particle # vanishes, as we have already seen earlier. 



B. Vorticity and superconductors 

As discussed in more detail in Paper I, the hydrodynamics of superfluids is generally characterized by two funda- 
mental properties: the absence of dissipative mechanisms like friction or viscosity, and the constraint of irrotational 
flow. While in the case of uncharged superfluids this simply meant the vanishing of the dynamical vorticity W, it 
is now the canonical vorticity W that is constrained to vanish identically in the case of charged superfluids, more 
commonly referred to as "superconductors". The absence of microscopic dissipative mechanisms implies that there is 
no non-conservative force acting on the bulk 2 of the superfluid, i.e. 

5 s = , (76) 

which quite generally characterizes perfect conductors of any sort. As a consequence we see that the canonical vorticity 
(and equivalently circulation) of a perfect conductor is strictly conserved, as seen in the previous section. The further 
constraint of irrotational flow, which distinguishes a superconductor from a mere perfect conductor, reads as 

w s = w s + -q s dA = 0, and W s = W s + -q s B = . (77) 

c c 

We see from ifTfljl or l(7T|l that if this irrotationality constraint is satisfied at some instant t, then it will automatically 
remain true for all subsequent times due to the absence of dissipation ifTfijl . We can therefore write the superfluid 
momentum 7r s (locally) as the gradient of a phase <p, i.e. 

a s 

tt s =p s + — A = HVu>, (78) 
c 

which leads to the well-known London equation for superconductors, as further discussed in Sect. IV Cl The canonical 
circulation l|73ll can therefore be non-zero if <9E encloses a topological defect in the phase <p, i.e. a region where (p 
(and therefore 7r s ) is not defined, as for example in the case of flow inside a torus, or around a vortex. While in the 
case of a perfect irrotational fluid the resulting circulation could have any value, the superfluid phase ip is restricted 
to change by a multiple of 2n when following a closed loop inside the superfluid around the defect. The resulting 
canonical circulation is therefore quantized as 

C = 2Nirh , with N e Z , (79) 

which gives rise to the well-known quantized vortex structure and flux quantization of superconductors. 



2 As mentioned in Paper I, this condition can be violated in the core of vortices, leading to "mutual friction". 
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C. Generalized helicity conservation 

We now turn to the generalization of the dynamical helicity conservation derived in Paper I. We define the canonical 
helicity 3-form H_ as the exterior product of the momentum 1-form n with the vorticity 2-form m, i.e. 

H = nAw, (80) 

and we define the dual canonical helicity density h as 

H = he, (81) 

where e is the volume form with components e^fc- The helicity scalar can be seen to have the following explicit 
expressions 

h = ttJW = 7T • (V x tt) . (82) 
Using l(fi9ll and (f7fil) , the comoving time-derivative of H_ can be found as 

{d t + £ v )H = d{Qm) + [d(z A g) + 2dl A tt] . (83) 
If we further introduce the total canonical helicity H of a volume V as 

U = I H= I hdV, (84) 
Jv Jv 

then we find in the absence of non-conservative forces, i.e. — 0, that the comoving time derivative of H satisfies 

= I QW-dS. (85) 

dt Jgy 

Therefore the canonical helicity Tt of a volume V is conserved under transport by the fluid only if, in addition to 
2 = 0, the canonical vorticity W vanishes on the surface dV surrounding this volume. We note that in general the 
conserved helicity H contains contributions from the purely hydrodynamic "Moffat" helicity p ■ W and the magnetic 
helicity A ■ B together with "mixed" terms, namely using 1821 we can express 

2 

h = p- W + %A-B + ^[p-B + A-W] . (86) 

cr c 



V. APPLICATIONS 



A. General description of electric conductors 

As a simple application of the foregoing formalism, we consider an electric conductor describable as a two-constituent 
system. One constituent consists of the positively charged ions, described by their number density n, velocity v, mass 
per ion m and charge per ion q = Ze. The second constituent is a gas of electrons of density n e , velocity v e , mass m e 
and charge q c = —e. The total charge density and current (pj are therefore expressible as 

a = e(Zn — n c ) , and j — e(Znv — n c v c ) , (87) 

and the relative velocity between the two fluids is 

A = v - v c . (88) 

Charge transfer between the two fluids is possible in principle, e.g. we could allow for processes of ionization and 
recombination, where electrons are transferred from the ion-fluid to the fluid of free electrons. But for simplicity we 
will assume the number of free electrons to be conserved, so we have 

r = d t n + V • (nv) = , and F e = d t n c + V • (n e v e ) = . (89) 

The total differential of the energy function £(n, n e , A 2 ) is 



d£ = /i dn + /i° dn c + a dA 2 



(90) 
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Using l-'t.'tt and ll-'Ut . the conjugate momenta of electrons and ions are therefore found as 

—A , 



p = mv 



-po = fi— \mv 2 + v ■ p , 



p c = m c v c + j^A , -pi = (j? - \m e vl + v c • p c 



(91) 



We neglect the gravitational field, so / gr£ 
are obtained from ll28l) and l(29jl as 



0, and the canonical force densities acting on the electron- and ion-fluid 



f = f a -nZe(E+-xB), and f c = f c H + en c {E- 



x B), 



(92) 



where the hydrodynamic force densities are obtained from ll26|l . by substituting the dynamical momenta (j9l"l) . which 
yields 



/ H = nm(d t + v ■ V) 
f% = n c m c (d t + v c -X7 



2a 



nV/i - 2aAjVv J 



2a 



n c m c 



2aAjVv J e . 



(93) 
(94) 



These equations contain the description of superconductors, magneto-hydrodynamic and the fluid-description of plas- 
mas (e.g. see @, El) as special cases. However, they are substantially more general due to the inclusion of the effect 
of entrainment, which is usually overlooked in these contexts. 

Using the momenta l(9"Tl) and the energy differential liffll the generalized pressure differential ijlfljl is found as 



c?vj/ = nd/i + n c d/i c — a dA 2 



(95) 



We note that in general we cannot introduce "partial pressures", say, by defining dP e to be equal n e d/j, e , as this is gener- 
ally not a total differential due to interaction energies between the constituents (i.e. the fact that fx e = /x e (n, n e , A 2 )). 
However, the chemical potentials are always well-defined and are therefore much more natural quantities in general 
multi-fluid contexts. In the absence of external forces, i.e. / cx t = 0, the force balance equation 14111 now reads as 



= / + / e = / H + /h - - J - x B . 



(96) 



We can further prescribe a mutual force between the two fluids, so we introduce a resistivity force of the form f c = fa, 
and therefore / = —Sr. The energy equation ijSOj) with l|24*j) now takes the form 



.9 + 9 C = - A • f R = g c 



(97) 



Such a resistive force will lead to creation of heat (entropy) , which in this model has to be extracted to an "external" 
system via g e xt, as for simplicity we have not explicitly included an entropy constituent in this example. By the 
second law of thermodynamics, the friction should produce heat and not absorb it, so we have to extract heat-energy 
from the system, i.e. g cx t < 0, and therefore we can constrain the resistivity force to be of the form 



f R = T]A, with ri > , 
where rj is generally a function of the state-variables describing the system. 



(98) 



B. The MHD limit 



In the low-frequency, long-wavelength limit we can assume any net charge densities to be compensated very quickly 
by the motion of electrons, so we make the "quasi-neutral" approximation and set: 

cr = , (99) 

which by Q87J1 implies Zn = n c , and the current density therefore reads as 



j = en c A 



(100) 
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In this low frequency limit we can equally neglect the displacement current dtD in Maxwell's equations 1181 . Because 
the electrons are very light, i.e. m c <C m, the inertial forces of the electron fluid can usually be neglected as well, and 
so the equation of motion for the electrons, f c — fn, can be written with plf . (IMI) and ijflSjl as 

n c Vfi c + 2aAjVvi + en c (E + — x B) = <qA . (101) 

In order to recover the "standard" MHD framework, we further neglect entrainment, i.e. if we set ad hoc a — 0, so 
the pressure differential ll95)) now reduces to 

dP = n c d^ c +ndfj,, (102) 

where the generalized pressure ^> can be identified with the usual pressure P in the absence of entrainment. This 
allows us to write the force balance equation (|96l) in the form 

p(9 t + ))-V)t) + VF--xB = 0. (103) 

c 



Using (flOOp we can express the electron velocity as 



1 , X 

v c = v j , 104 

en c 

and so we can write the equation of electron transport 11011) further as 

j = c(E + - x B) + -V// —j x B , (105) 

c e en e c 

where the scalar conductivity c is related to the resistivity coefficient r\ as 

c= ^£ >0 . (106) 
V 

If we further neglect the "partial pressure" V/i , we can write the relation between current j and electric field E' in 
the frame of the ion-background, i.e. E 1 = E + (v/c) x B &s & generalized Ohm's law, namely 

ji = c tk E k ' , (107) 

where the anisotropic conductivity tensor Cik is 

"1 „ 1 



which is not symmetric but satisfies the relation 



Sik H tikiB 1 

c en c c 



(108) 



c ik (B) = c fci (-B). (109) 

In this form the generalized Ohm's law can account for the well-known (classical) Hall effect, while the standard MHD 
approach (e.g. see ficjllll p commonly also neglects the "Hall term" j x B, so that this equation finally reduces to the 
standard Ohm's law: 

j = cE' = c(E + -v x B) . (110) 

c 

We note that the "orthodox" equations are contained in this framework as special cases, but the description (flDip 
is substantially more general due to the inclusion of the entrainment effect, which will generally be present in any 
(interacting) multi-fluid system. 



C. Superconductors 



In contrast to the previous application, superconductors are perfect conductors, so the electrons can flow past 
the ions without friction, i.e. rj = in i|98ll . Therefore we cannot neglect the inertial and pressure forces of the 
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electrons a-priori. As mentioned previously (cf. Sect. lIV Bj) . in addition to the absence of friction, superfluids are also 
constrained to be irrotational, so 

^•=V [< 7^=0, (111) 
and in its dual formulation this explicitly reads as 

W° = Vx7r c = Vxp° — —B = , (112) 

c 

which we will see after translation to the "orthodox" language represents the (second) London equation. In the absence 
of "external" forces acting on the electron fluid, using 12311 we can reduce the equation of motion for the electrons, 
/° = 0, to the form 

= <9f7r e — V7Tq = dtp — Vpq + e(VAo — —dtA) , 

= d tP c ~ Vpl + eE, (113) 

where we have used ill oil fit and the definition 10 of the electric field E. This equation is the (first) London equation 
and expresses the acceleration of electrons under gradients of their "potential" Pq and an electric field. This equation 
also guarantees that the constraint H112J1 remains automatically satisfied under the evolution of the electron fluid. 

The two equations 111211 and 111311 were originally proposed (albeit in their "orthodox formulation") by F. and 
H. London [12] and have been very successful in describing the phenomenology of superconductors, and in particular 
the behavior in electric and magnetic fields. 

We conclude this section by a "translation" into the orthodox formalism (cf. the discussion of superfluid 4 He in 
Paper I) . Using the entrainment relation I|91J1 we can express the electron momentum p c as 

n c 1 

— = v j , (114) 

m c ens 

where we have introduced the orthodox pseudo-density ns of superconducting electrons, namely 

7is = z , with e = . (115) 

1 — £ m c n c 

With this relation, Eq. H112J1 can now be written in its conventional form as 

Tn 

B = -cV x (Xj) , with A = , (116) 

where we have assumed that the background of ions is stationary and non-rotating, so V x v = and dtv = 0. 
Therefore Eq. H113J1 can now be written as 

E = d t {kj) + \vpli (117) 

where the "partial pressure" term Vpp = — V(/j c — |m c » c 2 + v c ■ p c ) is often neglected. Eqs. illlfi|l and <|1 1 T|l represent 
the orthodox formulation of the classic London equations as usually found in the superconductivity literature (e.g. 

see HE!!]]!)- 
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APPENDIX A: APPROXIMATE GALILEAN-INVARIANCE OF ELECTRODYNAMICS 



As pointed out in the introduction, the combined electro-hydrodynamics constructed in this paper makes no claim 
at being strictly Galilean- invariant. The underlying framework should be thought of as a fully (locally) Lorentz- 
invariant theory (as developed in of which we consider only the small- velocity, low-frequency regime up to 

and including effects of order (c _1 ). It is well-known (e.g. see [Ty|) that the first post-Newtonian corrections to 
uncharged particle-mechanics in a gravitational field are of order (c~ 2 ) , and the corresponding limiting theory is the 
strictly Galilean-invariant classical mechanics. This is not the case for electrodynamics. As pointed out in [jj, one 
cannot obtain a Galilean-invariant limit and keep the full Maxwell-equations, except for introducing the infamous 
ether. However, as we will show here, by restricting ourselves to a suitable "non-relativistic" regime of small velocities 
and low frequencies, the full framework of electrodynamics admits an "approximately" Galilean invariant formulation 
up to and including order (c _1 ). This will still hold true for the combined electro-hydrodynamics, as there are no 
(c -1 ) effects entering from the mechanical sector of the theory. By consistently neglecting terms of order (c -2 ), the 
combined theory can be considered as "approximately" Galilean-invariant in this sense. 

The well-known Lagrangian density A EM of the electrodynamic field = 2Vr M A i ,i can be equivalently expressed 
in terms of the two (frame-dependent) vectors E and B, namely 

A EM = -^-F^ = ^{E 2 - B 2 ) , (Al) 

l07T 071" 

where E and B, defined for an observer w M , are given by 

Ei* = F^u u , = F Xl u v . (A2) 

Note that both vectors are purely spatial and orthogonal to u^. In the Lorentz-frame defined by they can therefore 
be identified with the common 3-vectors used in the "3+1" Newtonian language describing the electric and magnetic 
field. Conversely, the field-tensor F^ v is uniquely specified in terms of E and B, namely 

F = uAE+* (uAB) , (A3) 

where * denotes the Hodge duality operation with respect to the 4-dimensional Levi-Civita tensor e, and A stands for 
the exterior product. In the Lorentz-frame we can write this in components as 



p — 



( 


—E x 


-Ey 


-E z 


E x 





B z 


-By 


Ey 


-B z 





B x 


\E Z 


B v 


—B x 






(A4) 



Now consider the effect of the transformation to a frame K' moving with velocity V relative to the original frame 

1/2 

K. Introducing the symbols (3 = V/c, (3 — |/3|, n = f3//3, and 7 = (l — (3 2 ) , the corresponding Lorentz- 



transformation matrix A can be written as 



A°' = 7: 



A ', = A^" = -/3 7 n^ (A5) 
A*'j = (7 - l)n j n k + S jk . 

Using this, we can transform the frame-components of F M „ to K 1 and translate them back into E 1 and B' , which 
yields (e.g. see (lfl |l 

E\ = E {{ , E' ± = 7 (Ex + (3 x B) , (A6) 
B\ = B||, B' ± = 7 (B ± -/3xE) , (A7) 

where the parallel (y) and orthogonal projections refer to the boost-direction f3. 

The range of validity of the formalism developed in this paper is the "non-relativistic" small- velocity, low-frequency 
regime, so we introduce the small parameter e <C 1 and require that all characteristic velocities v satisfy \v\/c ~ (e), 
and that time-derivatives of field-quantities are small compared to spatial derivatives, i.e. dtB/c ~ dtE/c ~ (e). 
Consequently we need to restrict ourselves to small Lorentz-boosts, i.e. we assume (3 ~ (e), and so 

7=(l-/?T 1/2 = l + (£ 2 )- (A8) 
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We can see that the Lorentz-transformation l|A6IA7j) of the electromagnetic fields can now be written as 

E' = E + /3x B + (e 2 ) , (A9) 
B' = B -0 x E+ (e 2 ). (A10) 

As pointed out in [jj, these transformations do not form a strict Galilean invariance group, because they fail to be 
additive. But they do form an "approximate" invariance-group up to order (e 2 ) in the sense discussed above. Namely, 
combining two boosts, 1 and 2 , we find 

E" = E' + (3 2 xB' = E + {0 1 + f3 2 )xB + (e 2 ) , (All) 
B" = B'-0 2 xE' = B- (/9 X + 2 ) x E + (e 2 ) . (A12) 

Let us consider the effect of the transformation ijAQIAlDjl on the electrodynamics field-Lagrangian (|Alj) , for which we 
easily find 

A; M = -L (E' 2 - B' 2 ) = ^(E 2 - B 2 ) + (e 2 ) . (A13) 

In principle this would conclude our demonstration, as both the Maxwell-equations and the Lorentz-force law are 
derivable from this Lagrangian, but for completeness we will also discuss their explicit transformation properties. 

Note that in the conventional Newtonian "3+1" language, the boost V results in the following transformations 
(using the fact that a; = ct) up to order (e 2 ): 

t' = t, (A14) 
x' = x-Vt, (A15) 
v' = v-V, (A16) 
d' t = d t \ x ,=d t + V-V, (A17) 

which are just the usual Galilean boost transformations. Similarly, we obtain the transformation law up to (e 2 ) for 
the gauge-field vector = (Aq, A) as 

A ' = A Q +(3-A, (A18) 
A' = A + Aa/3, (A19) 

and one can easily verify using the definition © of E and B in terms of (A , A) that this is consistent with the 
transformations ilA9IA10|) up to corrections of (e 2 ). It is also interesting to note the transformation properties of the 
charge-density a and the electrical current j. For an individual constituent X, the 4-current = nxu^, with the 
4- velocity — (c, vx) + (e 2 ). We therefore find up to order (e 2 ) 

°'x = °x , fx =3x-axV, (A20) 

and by the definition 10 of the total charge-density a and current j one easily finds 3 

a' = a , j'=j-aV. (A21) 

Applying the transformation ijAQHAlOjl and (IA 1411 to the Maxwell-equations l|17|l . I|18J1 . we obtain 

V • B' = p ■ (V x E) = -13 ■ (d t B/c) = (e 2 ) , (A22) 

VxE'l -diB' = V x E + -d t B -fix -d t E - (/3 ■ V)(/3 x E) + (e 2 ) 
c c c 



V xE+-d t B 

c 



+ (e 2 ) , (A23) 
V • E' - 47tct' = [V • E - 4ira] + V • (f3 x B) = -/3 ■ (V x B) = (e 2 ) , (A24) 



3 This shows that the current density always transforms according to the "electric limit" as defined in irrespective of the presence of 
net charges. The authors overlooked the fact that even a space-like current has to be the sum of time-like "elementary" currents, and so 
the "magnetic limit" will never apply in real systems. 
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I 47T I 1 

VxB - -d' t E' 3 = VxB - /3(V ■ E) - -d t E - (3 x -<9 t .B - (/3 • V)(/3 x B) 

c c c c 

47T 



+ (s 2 ) 



j - 4™/3 + (e 2 ) 



V x J3 - i<9 t £ - — j 

c c 



(A25) 



where in the last equation we used the transformation property of the electrical current density j' = j — aV . 
Finally, the expression for the Lorentz-force f28jl transforms as 



f 

J EI 



= E' + — xB' = E+-xB--x(3xE + (e 2 ) = — + (e 2 ) , 
nq c c c nq 

which concludes our demonstration of the approximate Galilean invariance of the equations. 



(A26) 
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